Weak-coupling approach to the semi-infinite Hubbard model: 
Non-locality of the self-energy 
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The Hubbard model on a semi-infinite three-dimensional lattice is considered to investigate 
electron-correlation effects at single-crystal surfaces. The standard second-order perturbation 
theory in the interaction U is used to calculate the electronic self-energy and the quasi-particle 
density of states (QDOS) in the bulk as well as in the vicinity of the surface. Within a real- 
space representation we fully account for the non-locality of the self-energy and examine 
the quality of the local approximation. Numerical results are presented and discussed for 
the three different low-index surfaces of the simple-cubic lattice. Compared with the bulk 
significant differences can be found for the top-layer local self-energy, the imaginary part 
of which is energetically narrowed and has a reduced total weight. The non-local parts of 
the self-energy E^(i£) decrease with increasing distance between the sites i and j. At the 
surface and for the three-dimensional bulk their decrease is faster than for a two-dimensional 
lattice. For all surfaces considered the effects of the non-local parts of the self-energy on 
the QDOS are found to be qualitatively the same as for the bulk: The weight of the quasi- 
particle resonance at the Fermi energy is lowered while the high-energy charge-excitation 
peaks become more pronounced. The main structures in the layer-dependent spectra are 
already recovered within the local approximation; taking into account the nearest-neighbor 
non-local parts turns out to be an excellent approximation. Due to the reduced coordination 
number for sites at the very surface, the top-layer QDOS is narrowed. Contrary to the the free 
(U = 0) system, quasi-particle damping results in a comparatively weak layer dependence of 
the QDOS generally. Pronounced surface effects that are related to surface Friedel oscillations 
show up as a fine structure within the resonance around the Fermi level. 

PACS: 71.27.+a, 73.20.At, 75.10.Lp 



I. INTRODUCTION 

The Hubbard model § g | is one of the most ex- 
tensively studied models for interacting fermions on a 
lattice. While originally it was set up to study the con- 
ditions for metal-insulator (Mott) transitions and spon- 
taneous magnetism, it nowadays quite generally con- 
tributes to the basic understanding of different types 
of prominent correlation effects. 

The mutual influence of electron correlations and the 
geometry of the underlying lattice is one of the central 
aspects of recent research. The simplest example for the 
decisive role of the geometry is given by the lattice di- 
mension D. In the extreme cases D = 1 (J, M, |(| and 
D = oo (j], H § the Hubbard model exhibits completely 
different physical properties. Also for a given dimension, 
say D — 3, the lattice structure is important. The exact 
results by Nagaoka H] as well as approximate treat- 
ments [O, represent illustrating examples for the 
interplay between magnetic order and the lattice type. 

The limit D = oo is of special interest since here the 
effects of the lattice geometry are suppressed as much 
as possible while the model itself remains non-trivial. 
The on-site Green function Gu a (E) depends on the ge- 
ometrical structure merely via the free Bloch-density 



of states (BDOS), and the self-energy T, a (E) becomes 
k-independent or site-diagonal (ll). However, even for 
D = oo the dependence on the lattice type, i. e. on the 
BDOS is important as has been proved, for example, 
by the recent investigations on ferromagnetism for the 
hypercubic and an fee- type lattice |lj], [T(| . 

It has been argued S that the essential physical prop- 
erties of the D — oo Hubbard model are comparable to 
those for finite dimensions. A solution of the D = oo 
Hubbard model would thus provide for a proper mean- 
field theory. Actually this argument can be considered 
as a justification of the local approximation for the self- 
energy, i. e. the neglect of its k-dependence or its non- 
local parts in the case of finite dimension. Since the 
contribution of non-local parts of the self-energy to rele- 
vant physical quantities scales with 1/D (|] , the local ap- 
proximation becomes less accurate in the opposite limit 
D i — > 1 . Starting from the local approximation and suc- 
cessively including nearest-neighbor, next-nearest neigh- 
bor and further correction terms, corresponds to an ex- 
pansion in 1/D. An instructive demonstration of this 
fact has been worked out by Schweitzer and Czycholl 
fTij] for D = 1 — 3. The expansion converges even for 
D = 1 ; for D = 2 all corrections beyond the 5-th neigh- 
bor shell are found to be negligibly small, and for D = 3 



the local approximation is quite appropriate. 

The interplay between geometric and electronic struc- 
ture becomes more important and interesting when the 
translational symmetry of the lattice is incomplete. This 
situation is realized e. g. in the case of a film of finite 
thickness composed of two-dimensional layers or at the 
surface of a three-dimensional lattice. Our study aims 
at the latter case, the semi-infinite Hubbard model. The 
semi-infinite Hubbard model may be interesting for dif- 
ferent reasons. First of all, it comes closer to the exper- 
imental situation. Because of the finite inelastic mean 
free path of low-energy electrons in the solid Q , sur- 
face effects are often non-negligible or even dominating 
in electron spectroscopy. Furthermore, surface effects 
may be interesting of their own. Particularly, the semi- 
infinite Hubbard model can contribute to the under- 
standing of phenomena related to surface magnetism. 
The combined influence of electron-correlation as well 
as of surface effects is important for the study of mag- 
netic surface states, of magnetic moments at the surface 
and a possibly enhanced surface Curie temperature etc. 
To our knowledge little effort has been spent on these 
topics up to now; the investigation of itinerant elec- 
tron models with reduced translational symmetry (thin 
films, surfaces) is still at the very beginning (see Refs 
" l2~ 
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23]), On the other hand, there has been 
extensive work on the semi-infinite Ising and Heisenberg 
model |24|]. These localized spin models, however, have 
to be considered as effective models where the magnetic 
order is more or less pretexted rather than derived con- 
sistently from the electronic structure. 

Despite its conceptual simplicity, an exact solution of 
the Hubbard model or an approximation scheme that 
is reliable in the entire parameter space is not avail- 
able for D = 3 up to now. The presence of the surface 
and the corresponding breakdown of translational sym- 
metry in the surface normal direction still complicates 
the problem, and thus approximations must be toler- 
ated. One of the main purposes of the present paper 
is to contribute to a clarification of the status of the 
local approximation for the self-energy in the context 
of the semi-infinite Hubbard model. This promises to 
be an interesting question since a priori it is not ob- 
vious whether a solid surface in this respect should be 
regarded as being close to D = 2 or close to D = 3. 
However, while the local approximation may be suffi- 
cient in the case D = 3, it is less accurate in D = 2 
fl7| , HI |6| . It can be expected that the degree of non- 
locality of the self-energy will be tightly connected with 
the reduction of the coordination number of sites at the 
very surface. One should try to substantiate a reason- 
ing based on coordination numbers by means of a more 
quantitative analysis. This includes the investigation of 
different (low-index) surfaces of the same lattice. While 
non-local parts Y,ij a (E) for i ^ j are expected to de- 



crease for distances R, — R, 



oo, the number of 



atoms within a shell of increasing radius R around the 
site i may increase. This raises the question of the effect 



of non-local contributions on the quasi-particle density 
of states and the changes near the surface. 

In the present paper we try to tackle the mentioned 
problems within the weak-coupling approach to the 
semi-infinite Hubbard model, i. e. within second-order 
perturbation theory (SOPT) in the interaction U. For 
the translationally symmetric Hubbard model this is a 
standard approximation scheme jlT], 27, ^8]. 29 31, 
[32[ |33| , p% [3f|. Its numerical evaluation in k-space in- 
cludes non-trivial multiple integrals. An efficient real- 
space algorithm for the calculation of the self-energy 
Y,ij a (E) was introduced in Refs. jl7|, [36] and is easily 
seen to be applicable also in the case of reduced trans- 
lational symmetry. The SOPT will yield exact results 
for infinitesimal U = + and will give first hints at least 
for finite U. An extrapolation of the results to strong 
interaction U (compared with the Bloch-band width W) 
must be questioned, of course. On the other hand, the 
comparison of the SOPT results for moderate interac- 
tion U [ |33| w ith the results of a strong-coupling moment 
method |lf[ is encouraging since qualitatively both ap- 
proaches agree well. In Ref. |33) Bulk and Jelitto em- 
ployed the SOPT around the Hartree-Fock (HF) solu- 
tion. This approach turns out to be superior compared 
with the plain SOPT in which free instead of Hartree- 
Fock propagators are used in the calculation of the 
proper irreducible self-energy diagram. Neither plain 
SOPT nor SOPT-HF, however, is a conserving approx- 
imation in the sense of Baym and Kadanoff |37j . Con- 
trary, the self-consistent SOPT (using the fully renor- 
malizcd propagators) does represent a ^-derivable the- 
ory H g. 

Most of the results presented here have been calcu- 
lated for U — + . In this case the different versions of 
the weak-coupling theory become identical. Consider- 
ing an infinitesimal interaction, however, does not allow 
to study the manifestation of correlation effects in the 
quasi-particle density of states. Therefore, we will also 
consider finite interaction U, but restrict ourselves to 
the symmetric situation of a paramagnet at half-filling 
(n = 2 (71^) = 1) in this case. The symmetric case 
is somehow exceptional; the following well-known ar- 
guments substantiate that SOPT-HF may yield qual- 
itatively correct results even for moderate interaction 
U ~ W at n = 1. Firstly, SOPT-HF exactly reproduces 
the atomic limit (W = 0) at half-filling. Secondly, the 
energetic positions as well as the spectral weights of both 
high-energy charge-excitation peaks that show up in the 
spectrum for U ^> W can be calculated exactly within 
a strong-coupling perturbation approach (expansion in 
W/U) m |§. For the case of half-filling SOPT-HF can 
easily be shown to reproduce these exact results. We 
thus believe that the method reasonably interpolates 
between the weak- (U <C W) and the strong-coupling 
(U 3> W and W = 0) regimes. On the contrary, the 
self-consistent SOPT fails with respect to all points men- 
tioned. In our opinion this cannot be outweighted by 
its advantages, namely to represent a conserving theory 



and thus to respect the Luttinger sum rule [}43| . Since 
the SOPT-HF scheme preserves particle-hole symmetry 
|Q , the sum rule is fulfilled anyway in the case of half- 
filling. The last argument regards the high-energy be- 
havior of the SOPT-HF self-energy. At half-filling the 
first two coefficients (for the local part three coefficients) 
within an l/E expansion turn out to be exact (see Ref . 
]4l| , ffi|). This also implies that the first three moments 
of the quasi-particle density of states can be reproduced. 
Another reason for using SOPT-HF is a pragmatic one: 
For the general (non-symmetric) case, self-consistency 
with respect to the occupation number (n icr ) and the 
chemical potential p is required. At half-filling this is 
achieved trivially, we have p — U/2 (this also holds for 
the semi-infinite lattice). Since our main interest is fo- 
cussed on the validity of the local approximation for the 
self-energy, we may set aside the fact that the Hubbard 
model at half-filling is expected to be an antiferromag- 
net Q. Throughout the paper we exclusively consider 
the (possibly metastable) paramagnetic phase. 

The next section briefly introduces the theoretical 
concept. In section III we present and discuss the re- 
sults. All calculations have been performed for zero tem- 
perature T = 0. We have limited our investigations to 
the simple-cubic (sc) lattice. In most cases the (100) 
surface was considered, but also the more open sc(110) 
and sc(lll) surfaces were taken into account for com- 
parison. A summary and the conclusions are given in 
section IV. 

II. THEORY AND CALCULATIONS 

We consider interacting electrons moving through a 
semi-infinite lattice. The dynamics of the system is de- 
scribed by the Hubbard model. Using standard nota- 
tions the Hamiltonian reads: 



(1) 



The system geometry is implicit: c* (c) is the creation 
(annihilation) operator for a valence electron. The lower 
indices specify the orbital of the one-particle basis, a =f 
, J, accounts for the spin direction, and i and j refer to 
the sites within the semi-infinite lattice. We consider 
the system to be built up from single two-dimensional 
layers parallel to the surface. Accordingly, the site index 
may be decomposed as i = (i±,i\\) where i± labels the 
layers i± = 1, . . . , oo and i\\ the sites within a given layer. 
i± = 1 stands for the topmost surface layer. Within 
a tight-binding approach the hopping integrals are 
assumed to be non-zero up to nearest neighbors only 
and constant throughout the system, i. e.: 2y = — t 
for nearest neighbors i, j. The energy scale is fixed by 
taking 4 = 1. The energy zero will be chosen such that 
Tu = To = 0. U denotes the on-site Coulomb interaction 
which is assumed to be layer independent. A change of 
U as well as of T^- in the vicinity of the surface should 



be expected for any real system |4fi(] . Here the constancy 
of the parameters is a model assumption which helps to 
focus on the most important geometrical aspect, namely 
the reduction of the coordination number at the surface. 
Finally, in Eq. ([j]) n,- CT = c^Cfo is the local occupation- 
number operator, and p denotes the chemical potential. 

All properties of the system that we are interested in 
can be obtained from the one-electron Green function: 



G ija {E) = ({c irT ;c\ )) J 



(2) 



Let us assume the one-particle basis orbitals to be real. 
We then get the usual diagonal i = j as well as the off- 
diagonal i j quasi-particle density of states (QDOS) 



via: 



Pi jcT {E) = lmG ija (E + iO^ 



(3) 



For the semi-infinite system the QDOS will be layer de- 
pendent in the vicinity of the surface. The Green func- 
tion obeys an equation of motion which in real-space 
representation reads: 



EG ija (E) = hS %3 +J2 (?ik - pfa + ^ika(E)) G kja {E) . 

k 

(4) 

Here we have introduced the electronic self-energy £ 
which includes all effects of electron correlations. 

Within the usual perturbative approach p7| the self- 
energy is expanded in powers of the interaction U . The 
first non-trivial term is proportional to U 2 . Using the 
real-space notation again, we have: 



Z ija (E) = U(n i - <T )5 ij + U 2 EgJ(25) 



(5) 



The term linear in U is the Hartree-Fock self-energy. 
Restricting ourselves to zero temperature, the second- 
order contribution reads: 



1 



E — x + y — z 
(&(-x)e(y)Q(-z) + @(x)G(-y)G(z)) dx dy dz . 

(6) 

8 denotes the step function. Within plain SOPT is 
the free (off-diagonal) density of states: 



p$l(E) = -~ImG ijlT (E + iO+: 



(7) 



The SOPT-HF approach (cf. e. g. Refs. 0, |§) replaces 
the free by the Hartree-Fock density of states in the 
expression (||): 



lmG?){E + m+), 



(8) 



where the HF Green function G^ has to be calculated 
from 



EG%l(E) = h5 ij 

+ E { T '" + ( U ( n i-°) - M) &*) G kl(E) . (9) 

k 

The ground-state expectation value of the local particle 
number is obtained by integrating the (fully interacting) 
QDOS: 



1 f° 

(»»(t) = ^ J Piia{E) dE . 



(10) 



Eqs. (||) - ( ^0| ) constitute a closed set of equations 
for the second-order perturbation theory around the 
Hartree-Fock solution. Generally, the local particle 
numbers (n icr ) have to be determined self-consistently. 
Due to the presence of the surface, the self-consistent 
values are expected to be layer dependent, i. e. depen- 
dent on the distance from the surface. A self-consistent 
solution is not necessary in two cases, namely for an in- 
finitesimally small interaction U — + and for the sym- 
metric case of a paramagnet at half-filling. In the latter, 
we have {riia) — (ni_ ff ) = 0.5 for all particle numbers, 
and the chemical potential is given by /i = U/2. (In 
the following the spin index will be omitted as far as 
possible; throughout the paper only paramagnetism is 
considered) . 

There are two main problems associated with the nu- 
merical solution of the SOPT-HF equations. The first 
one consists in the non-locality of the second-order term 
of the self-energy (^j). For the calculation of the local 
contribution i = j a three-dimensional energy integral 
has to be performed. Considering the imaginary part 
of the self-energy Im £^?' (E + i0 + ), this reduces to con- 
volution integrals which can be calculated efficiently by 
applying a Laplace transform. The real part is then ob- 
tained through the Kramers-Kronig relation. Details of 
this procedure can be found e. g. in Refs. [jl^, ^|. Ne- 
glecting all off-diagonal elements in (0) corresponds to 
the local approximation for the self-energy. The numeri- 
cal effort needed to account for the complete non-locality 
of the self-energy can be estimated in the following way: 
Let us consider the neighbor shells around the site i. If 
there was perfect translational invariance with respect 
to all three spatial dimensions, symmetry reasons would 
require all non-local terms to be equal for sites j be- 
longing to the same shell around the site i. Therefore, 
compared with the local approximation the effort rises 
by a factor s + 1 that is determined by the number of 
shells s that are necessary to obtain convergence for the 
results JT^j (as mentioned in the introduction, the neces- 
sary number of shells strongly depends on the lattice di- 
mension D). In the semi-infinite system the symmetry is 
reduced, i. e. not all the sites j within the same shell are 
equivalent. Consider, for example, a (fixed) site i near 
the surface. Independent from the lattice type and the 
surface there are three incquivalcnt possibilities at least 
to choose a site j within the nearest-neighbor shell of i: 



j within the same layer as i (J±_ — i±), j within a layer 
more distant from the surface = i±_ + n with n > 0), 
and finally j closer to the surface (j± = i± — n). (Note 
that nearest neighbors need not necessarily belong to 
adjacent layers.) Therefore, the reduced symmetry im- 
plies an increase in the numerical effort by a factor three 
or even more, since the number of incquivalent sites j 
within the same shell around i increases with increasing 
shell radius. Finally, we have to consider all inequivalent 
positions for the site i. Assuming perfect translational 
symmetry within the layers parallel to the surface, their 
number is given by the number of inequivalent layers 
N± . The computational effort therefore rises once more 
by a factor Nj_. With increasing distance of i and j 
to the surface, (E) converges to the bulk self-energy. 
Thus N± can be regarded to be finite but may be rather 
large (several tens of layers) . This also implies that in- 
stead of the semi-infinite we can assume a slab geometry 
with a large number of layers d (~ 2Aj_) in the actual 
calculation. 

The second main problem concerns the equation of 
motion ([|) for the Green function. It is the reduced 
translational symmetry of the semi- infinite system again 
that renders its solution more difficult. In order to 
exploit translational symmetry within the layers, two- 
dimensional Fourier transformation of any quantity Fij 
will be convenient: 



F i4 = F 



) F^^A\\), (11) 



and 



Fs (kill = — V e _ik|1 (R, ii " Rj ii } F' i±3j 
11 *uiii 



(12) 



Here k|| is a wave vector of the first two-dimensional 
surface Brillouin zone. After Fourier transformation the 
equation of motion for the Green function reads: 



EG ixj± (k h E) = fiS i±j± +J2( T ixkAh) 



-fJ^ixkx + ^ixkx ( k || ' E )j G kx3X ( k || ) E ) ■ (I 3 ) 

This may be considered as a matrix equation for Gi ± j ± . 
The dimension of the inversion problem is given by the 
slab thickness d ~ 2N±. Within the tight-binding ap- 
proach, the hopping matrix Ti ± j ± as well as the self- 
energy matrix ^i x j ± are sparse. Therefore, the tight- 
binding recursion technique 48 , 49 2l[] is most suitable 
for its solution. The application of the method yields 
the Green function as a continued fraction: 



Gi ± i ± (k-\\, E) — G l±l± — G[ ^ , 



G 



(fc) 



<fc) _ ( h {k+i)Y G (k+i) 



e -<:-[Kx 



(14) 



The (energy-dependent) coefficients a\ k ' and b^ k+1 ^ can 
be calculated iteratively within the recursion approach. 
Due to the finite dimensionality, the continued fraction 
terminates automatically: k = 0, . . . d — 1. Off-diagonal 
elements of the Green function can be obtained by means 
of a polarization relation [|TJ . 

III. RESULTS 

We have chosen the simple-cubic (sc) lattice as a 
model system for our calculations. Three different low- 
index surfaces, sc(100), sc(110) and sc(lll), have been 
considered. Let us start the discussion of the results with 
the (100) face. For the sc lattice the nearest-neighbor co- 
ordination number for sites in the bulk is z — 6. Speak- 
ing in terms of the (100) layers parallel to the surface, 
each site has z^ 100 ^ — 4 nearest neighbors within the 

same layer and z^ 100 ^ = 2 sites in the adjacent layers. 
For the top-layer sites (i± = 1) the coordination num- 
ber is thus reduced to Zg 100 ' = 5. 

This (small) lessening of the coordination has conse- 
quences for the free (U — 0) density of states (DOS) 
pf^{E). The free DOS determines via Eq. || the second- 
order contribution of the self-energy for U = + . For fi- 
nite U and at half-filling we have to consider the Hartree- 
Fock (HF) DOS p$(E). Compared with the free DOS 
the HF-DOS is only shifted constantly in energy by the 
amount U(rii- a ) = U/2 which is compensated by the 
corresponding shift of the chemical potential. So in both 
cases pfj' ] (E) is the primary quantity of interest. 

Let us discuss first the layer-dependent diagonal el- 
ements (i = j) of the free DOS which are shown in 
Fig. 1. We notice that there are considerable differences 
between the DOS of the topmost surface layer and the 
bulk DOS. The difference becomes weaker and weaker 
with increasing distance from the surface. The surface 
can be regarded as a perturbation of the infinitely ex- 
tended lattice that gives rise to oscillations of the layer- 
dependent charge density. These surface Friedel oscilla- 
tions J50) are closely related to the periodic deformations 
of the DOS which, as can be seen in Fig. 1, are damped 
when passing from the very surface to the bulk. The 
layer-dependent off-diagonal DOS pfj\E) for nearest- 
neighbors i and j within the same layer is also shown 
in Fig. 1. Its absolute magnitude is smaller, and there 
is a zero of the DOS at the center of gravity Tq = 
of the diagonal DOS. Compared with the case i = j, 
we observe qualitatively similar differences between the 
top-layer DOS and the bulk DOS as well as the signature 
of surface Friedel oscillations. 

The overall shape of the free DOS in the vicinity of 
the surface can be understood in terms of the moments 
of the DOS. These are defined as: 

/oo 
E n p[f(E)dE. (15) 
-00 

An alternative, but equivalent representation can be de- 
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Fig. 1. Free (U = 0) densities of states p u for the first three 
layers from the surface (i± = 1 — 3) and free bulk density of 
states (i± = 00) as functions of energy (solid lines). i± = 1 
stands for the outermost surface layer. Dashed lines: off- 
diagonal free densities of states pf) (E) for nearest-neighbor 
sites i and j, both within the same (indicated) layer. Calcu- 
lations for the sc(100) surface. (A slab thickness of d = 50 
turns out to be sufficient for convergence in all calculations) . 
All energies are given in terms of the nearest-neighbor hop- 
ping t = 1. p = 0. 

rived from the equation of motion for the free Green 
function: 

M K0) = £ T Ul T ili2 ---T in _ d . (16) 

Using this equation we can conclude that for each layer 
the diagonal (i = j) DOS must be a symmetric function 
of energy with respect to its center of gravity Tq = 0: 
All odd moments must vanish since there is no way to 
start from and to return to a site i by an odd number 
of nearest-neighbor hoppings for a (semi-infinite) simple 
cubic lattice. Analogously, the off-diagonal free DOS for 
nearest-neighbors i and j must be antisymmetric. 

Comparing with the bulk, the on-site top-layer DOS 
is much more compressed. This "band-narrowing" is a 
consequence of the reduced coordination number of the 
top-layer sites and can be understood by considering the 
second moment: 

av1 0) = mT } - (M^ o) r = E T l = *«) t2 ■ ( l7 ) 

The reduced coordination number z(i) of the site i at 
the surface yields a reduced variance A 2 p(°*> of the top- 
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Fig. 2. Free off-diagonal densities of states p^(E) near the 
sc(100) surface. The respective sites i and j are indicated 
as filled circles in the corresponding schematic drawing of 
the surface geometry. The shaded planes represent the first 
three layers from the surface, i± = 1 — 3 (i± = 1: top surface 



layer). Note that p ( ° ] '(E) 



0. 



layer DOS. For the second and the third surface layer 
higher moments have to be considered. Despite the 
band- narrowing effect, the densities of states of all lay- 
ers share common band edges. For a spatially constant 
hopping parameter t no split-off states can be observed. 

The difference between the free densities of states for 
non-equivalent pairs of sites i and j is demonstrated in 
Fig. 2. For the calculations the site i has been kept 
fixed within the second layer from the surface (i± = 2). 
The first three panels from the top in Fig. 2 show the 
off-diagonal DOS pfj\E) for nearest-neighbor sites j. 
There are three non-equivalent positions relative to i. 
We recognize significant though not strong differences 
between the spectra. These are exclusively due to the 
presence of the surface; choosing for i a site within the 
bulk would yield identical results in all three cases. The 
last three panels show the DOS for i within the second 
layer as before, but j a next-nearest neighbor. Again, 
three non-equivalent positions for j can be found. Since 
an even number of nearest-neighbor hoppings is needed 
to get from site i to site j, the DOS is symmetric. The 
differences between the spectra are non-negligible and of 
the same order of magnitude as in the nearest-neighbor 
case. They will lead to a (weak) directional dependence 
of the self-energy in the vicinity of the surface. 

The local approximation for the self-energy Ey(i£) 
can only be a reasonable starting point if the non-local 
parts strongly decrease with increasing distance between 
the sites i and j. Within the perturbational approach 
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Fig. 3. Free densities of states p^'(E) near the sc(100) sur- 
face (solid lines). Sites i and j are indicated (shaded planes: 
(100) layers, top plane: surface layer). Dashed lines: bulk 
densities of states. Relative position of sites i and j as indi- 
cated (shaded planes now to be interpreted as (100) layers 
in the bulk), p = 0. 



these are determined by the off-diagonal free DOS. In 
Fig. 3 the DOS pf^ (E) is shown for a site i fixed within 
the second layer from the surface and different sites j 
with an increasing distance dij = ||Rj — Rj|| (from the 
top). It can be seen that the absolute magnitude of the 
DOS decreases with increasing distance dij indeed. Par- 
ticularly, the off-diagonal DOS for nearest neighbors is 
clearly smaller in absolute values compared with the on- 
site DOS over the whole energy range. A more quanti- 
tative estimate can be given by integrating the absolute 
DOS. We then obtain the ratio (for increasing distance, 
on-site DOS normalized to 100) 100 : 33 : 21 : 28 : 14 : 9. 
The decrease is not monotonous. According to Eq. 
(H) a stronger decrease has to be expected for the U 2 - 
contribution to the self-energy since £^ (roughly) re- 
sults from the third power of the DOS (p\^) 3 - 

The results for the site i located near the surface can 
be compared with the results for a site i in the bulk 
which are also shown in Fig. 3. The differences between 
the corresponding surface and bulk off-diagonal densi- 
ties of states are significant though not strong. They 
are qualitatively similar to the difference between the 
surface and bulk on-site DOS. Since the second-order 
contribution to the self-energy (0) results from a (two- 
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Fig. 4. Real and imaginary part of the infinitesimal on-site 
self-energy | -^yzY,u{E + iO + )|[7=o of the semi- infinite sc(100) 
Hubbard model. Band-filling: n = 0.6. Solid line: site i 
within the topmost surface layer (i± = 1). Dashed line: 
bulk self- energy. 
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(18) 

The difference is found to be small: w = 0.21 and w s = 
0.20. 

Both, the imaginary part of the top-layer and of the 
bulk self-energy, vanish quadratically, — ImS^(S) ~ 
E 2 , as E i — * 0. This follows from a straightforward 
analysis of Eq. (||) |il| and implies the existence of a 
well-defined Fermi surface fl52| . We can thus distin- 
guish between the damping of the occupied (E < 0) 
and of the unoccupied [E > 0) part in the excita- 
tion spectrum. Let and m?\ denote the integrated 

weight of — ~lmS!?\E + i0 + ) at the surface or in the 
bulk for E < and E > 0, respectively. We have 
w (s) w (s) = '"(b)) an d from Eq. (|g) the following re- 
sult can be derived: 



fold) convolution of the free DOS, these differences be- 
come more or less meaningless. We thus can conclude 

( 21 

that the self-energy Y^-' {E) for a site i within second 
layer and for an arbitrary site j is almost identical to 
the corresponding diagonal or off-diagonal self-energy in 
the bulk. Therefore, for a site i in the second surface 
layer the local approximation is likewise appropriate (or 
inappropriate) as for a site i in the bulk. 

Let us consider now the top layer where the situation 
is found to be different. Fig. 4 shows numerical results 
for the infinitesimal on-site self-energy i^j T E^(i?)|[/=o 
of the semi-infinite sc(100) Hubbard model. We con- 
sider a non-symmetric case: The chemical potential has 
been chosen such that in the bulk the average occupa- 
tion number (band-filling) is n = 2{ni a ) — 0.6. It can be 
seen in Fig. 4 that there are significant differences be- 
tween the top-layer and the bulk self-energy. Contrary, 
the infinitesimal on-site self-energy of the second layer is 
almost indistinguishable from the bulk self-energy. On 
the energy scale used in the figure, differences would 
not be visible. The same holds true for the third and all 
other surface layers. 

Fig. 4 shows that the on-site self-energy of the top 
layer is more compressed compared with the bulk self- 
energy; the effective width of the imaginary part is 
smaller at the surface. This is a typical surface effect 
since it is a direct consequence of the narrowing of the 
top-layer DOS mentioned above. Another consequence 
concerns the local occupation number (n irT ) at the sur- 
face. Due to the narrowing of the top-layer DOS, the 
surface charge density n s = 2{rii a } (for i within the first 
layer) must be smaller than the bulk band-filling: n s < n 
(below half-filling) @. For n = 0.60 we find n s = 0.56. 
This leads to a reduced total weight w s of the imaginary 
part of the infinitesimal self-energy at the surface com- 
pared with its weight in the bulk w: From Eq. (^3j) we 
have: 




It can be seen in Fig. 4 that for the occupied part the 
integrated weight is clearly smaller in the top layer. This 
has to be interpreted as another typical surface effect: 
The reduction of the coordination number at the surface 
implies the narrowing of the top-layer DOS and thus 
the lowered charge density; via Eq. (|ll]) this leads to a 
smaller damping in the occupied part. We get wf — 

0. 056 to be compared with the bulk value w < — 0.063. 
On the contrary, the weight within the unoccupied part 
is almost unchanged: w> — 0.145 and w > = 0.147. 

The non-local parts of the infinitesimal self-energy for 
the top layer of the sc(100) surface and for the bulk 
are shown in Fig. 5. Unlike the imaginary part of on- 
site (retarded) self-energy, which must be negative over 
the whole energy range, there is no definite sign for the 
imaginary part in the off-diagonal case. Instead, we have 
the sum rule /^ImS^B + iQ+) dE = for i ^ j. 
The low-energy behavior, however, remains unchanged: 
the imaginary part vanishes quadratically for E 0, 

1. e. ImS[^(£) ~ ±E 2 for i =/= j. The sign depends on 
the sign of the (off-diagonal) DOS at the Fermi edge. 
This can be verified by expanding the infinitesimal self- 
energy in powers of E. Similar to the analysis in Ref. 
Jslll , one may derive the following result in the real-space 
representation: 

-ilm£g>(£+z0+) = (jrpi 0) (0)) xE 2 +0(E% (20) 

Consider, for example, nearest-neighbor sites i and j and 
n < 1. As can be seen from Fig. 1, pf^{E) is positive 
at the Fermi energy and thus lmY$ (E + i0 + ) ~ — E 2 
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Fig. 5. Off-diagonal elements of the imaginary part of 
the infinitesimal self-energy i -^ I T,ij(E + i0 + ) \u=o near the 
sc(100) surface (solid lines) and in the bulk (dashed lines). 
Band-filling n = 0.6. Sites i / j as indicated (dashed 
lines/bulk: shaded planes represent (100) bulk layers). The 
vertical energy scale is given by the height of the bar (0.001 
in units of £). In the last two panels the data have been 
multiplied by a factor 10. 
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Fig. 6. Imaginary part of the (infinitesimal) self-energy 
i^ T S ij (£+i0+)|(7=o for the two-dimensional (D = 2) Hub- 
bard model on the square lattice at half-filling (n — 1). Site 
j within the fc-th neighbor shell of the (fixed) site i. Results 
for k = (on-site self-energy) up to k = 7. The results for 
k ^ have been scaled as indicated. 



for E i ► which is consistent with the result shown in 
Fig. 5. 

The first two panels from the top in Fig. 5 show the 
imaginary part of the self-energy for two non-equivalent 
pairs of nearest-neighbor sites i and j. A significantly 
stronger variation of the self-energy is observed when 
both sites are located within the top layer. This implies 
a directional dependence of the self-energy which, how- 
ever, is confined to the very surface. In the bulk both 
cases are equivalent, and the same result for the self- 
energy is found. Generally we can state that significant 
though not strong differences to the bulk self-energy are 
only found in those cases where both sites i and j belong 
to the topmost surface layer. This again proves that the 
presence of the surface manifests itself in the self-energy 
for the top surface layer only. 

Increasing the distance dij between the sites, results 
in a rapid decrease of the absolute values for the self- 
energy in the bulk as well as in the vicinity of the surface. 
Compared with the on-site self-energy (Fig. 4), the ab- 
solute values of the non-local self-energy are smaller by 
more than one order of magnitude already for nearest- 



neighbor sites i and j. We have integrated the energy- 
dependent expression |Im X \-' (E+i0 + )\ for i = j (Fig. 4) 
and for i ^ j (second to sixth panel in Fig. 5) to get a 
more quantitative estimate for the magnitude of the dif- 
ferent contributions at the very surface. Normalizing the 
local contribution to 100, the following ratio is obtained 
for increasing distance: 100 : 4.3 : 1.4 : 1.6 : 0.09 : 0.20. 
As had to be expected, this roughly corresponds to the 
above-mentioned ratio for the third powers of the inte- 
grated (absolute) densities of states. It shows up that 
the local part of the self-energy (bulk as well as surface) 
is rather dominating. 

The results for the semi-infinite lattice may be com- 
pared with those for a two-dimensional case. For this 
purpose we have calculated the local and the non-local 
parts of the infinitesimal self-energy on the D ~ 2 square 
lattice. Fig. 6 shows the results for the imaginary part at 

(2) 

half-filling. The self-energy S^ . (E) is a symmetric (an- 
tisymmetric) function of energy with respect to E = 
if an even (odd) number of nearest-neighbor hoppings is 
necessary to get from site i to site j. At half-filling the 
low-energy behavior of the self-energy is likewise excep- 



JD 

lET 

Cl 
CO 

3 
c5 

-t— » 
CO 

H— 

O 
>* 

4—* 

CO 
CD 

Q 



1 1 1 1 1 1 1 1 1 1 1 1 1 

U=8 


i i i i | i i i i | i i i 

D=2 


n=1 i 




c-D ^-^ J i 




S=1 I 




/ — I 




s=3 N^i 




/ 

A / 

s=4 IW 


\ 

\ h 

\ A 


s=5 A-^ 




i i i I i i i i I i i i i I i i i i I i i i i I i i i 



-10-5 5 10 
Energy E 

Fig. 7. Quasi-particle density of states p(E) of the D — 2 
Hubbard model on the square lattice at half-filling and U = 
8 within SOPT-HF. Off-diagonal contributions of the self- 
energy included up to the s-th neighbor shell. Results for 
s = (local approximation) up to s = 5. Dashed line: s — 
result for comparison. 



tional: Expanding the expression u 

1 ( 2) 

for the on-site part: — -Im£^ (-E 



for small E 1 yields 
i0+) - £ 2 |ln|£|| 3 
as B i— > |k| . The reason for this anomalous behavior 
at half-filling is the logarithmic divergence of the free 
DOS at E = in two dimensions. 

As can be seen from Fig. 6, the non-local parts of 
the self-energy are much smaller in absolute magnitude 
compared with the local one. Again they are found 
to decrease with increasing distance dij. For the in- 
tegrated absolute self-energies we obtain the following 
ratio: 100 : 10 : 7.4 : 1.1 : 2.1 : 2.0 : 0.80 : 0.52 (on- 
site part normalized to 100). Compared with the results 
found for the sc(100) surface, we notice a considerably 
slower decrease with increasing distance dij. Qualita- 
tively, this result does not depend on the band-filling at 
all. Although the mere numbers should not be overem- 
phasized, the comparison of the different ratios surely 
gives an impression about the relative importance of the 
non-locality of the self-energy in the different cases. 

So far only the free DOS as well as the self-energy 
have been considered. To be able to judge on the qual- 
ity of the local approximation, however, one also has 
to check the effects of the non-local parts on the quasi- 
particle density of states (QDOS): While |£$ ) (J5+iO+)| 
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Fig. 8. Layer-dependent quasi-particle density of states 
pu(E) for the semi-infinite sc(100) Hubbard model at half- 
filling and U = 8. SOPT-HF results for the first five layers 
from the surface and for the bulk. Dashed lines: local ap- 
proximation. Solid lines: non-locality of the self-energy fully 
taken into account. 



may be a small quantity for a given pair of sites i j, 
it nevertheless may be important with respect to the 
QDOS provided that a large number of equivalent pairs 
have to be taken into account. For the discussion of the 
(paramagnetic) QDOS at finite interaction U, we re- 
strict ourselves to the symmetric case of half-filling. As 
mentioned in the introduction, the SOPT-HF scheme 
reproduces a large number of exactly solvable limits for 
n = 1 and should thus yield reliable results, at least on 
a qualitative level. 

The QDOS of the D — 2 square lattice at half-filling 
is shown in Fig. 7. We have chosen U = 8 which is equal 
to the free Bloch band width W. The first panel shows 
the result for the local approximation (s = 0). The 
other panels show calculations where we have included 
the non-local parts of the self-energy up to s-th neigh- 
bor shell (s — 1 to s = 5). Qualitatively, all spectra look 
similar: The high-energy charge excitation peaks (Hub- 
bard bands) are clearly visible for U = 8, and a quasi- 
particle resonance remains around E = 0. Yet there are 
important differences: The inclusion of more and more 
shells narrows the whole spectrum, sharpens the charge- 
excitation peaks while the resonance looses some spec- 
tral weight. Furthermore, a new structure comes into 
existence at £ « il. Fully converged results cannot be 
obtained until s = 5 f54fl. 
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Fig. 9. Real part of the infinitesimal on-site self-energy 
^-^j I T,u(E + iO + )|iy=o for the top surface layer of the semi- 
infinite Hubbard model. Results for three different surfaces. 
Band-filling: n — 0.6. 



Let us now return to the sc(100) surface again. The 
layer-dependent QDOS in the vicinity of the surface and 
in the bulk for U — 8 and at half-filling are shown 
Fig. 8. For each layer the QDOS exhibits three main fea- 
tures which again can be interpreted as the two charge- 
excitation peaks (on the high-energetic sides) and a res- 
onance (at the Fermi level). 

The discussion of the self-energy above has already 
indicated that effects of its non-locality are less impor- 
tant for the semi- infinite system compared with D = 2. 
This becomes manifest when considering the QDOS at 
the surface. Taking into account the on-site self-energy 
and the nearest-neighbor non-local part is almost suffi- 
cient to obtain fully convergent results. Remaining dif- 
ferences between the QDOS for s = 1 and for higher s 
would hardly be visible on the scale used in Fig. 8. Let 
us point out that the convergence with respect to s is 
equally fast for both, the top-layer and the bulk QDOS 

The results shown in the figure prove that the local 
approximation provides a reasonable description of all 
features appearing in the spectra. Only the spectral 
weight of the resonance at the Fermi edge is somewhat 
overestimated, and the charge-excitation peaks are less 
pronounced for s = 0. Compared with the results for the 
D = 2 square lattice, however, the local approximation 
is only slightly better. In any case it is likewise appro- 
priate for the bulk as well as for the surface including 
the top layer. 

The results for the QDOS at the sc(100) surface are in- 
teresting of their own: First, we notice that the top-layer 
QDOS is slightly narrowed. This can be understood by 
referring to the moments again. The first few moments 
of the QDOS, 
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/oo 
E n p l3 (E)dE, (21) 
-OO 

can be calculated exactly (see Rcf. J2TJ , for example): 
Mi° } = 1 



M 
M 



(i) 



= T + U(ni^ a ) 



(2) _ 



J2T?j + 2UT (ni- a ) + U 2 (n^ a ) . (22) 



The variance of the QDOS for the layer i± is then given 
by the expression: 



A 2 p u = Mj?-(MXr = 



(1) ^ 2 -E^- 



-U 2 (n,, 



r )(l-(rw)). 

(23) 

We can conclude that increasing U leads to a strong 
increase of the variance of the QDOS which is mainly 
due to the evolution of the two Hubbard bands. For a 
given U and at half-filling the second term on the right- 
hand side of ( p3[ ) is a constant. As for the case of the 
free DOS, the first term implies a narrowing of the top- 
layer QDOS because of the reduced coordination number 
for sites in the first layer. This can be seen in Fig. 8. 
Additionally, we notice a slight narrowing of the quasi- 
particle resonance for the top layer. 
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Fig. 11. Layer-dependent quasi-particle density of states for 
the sc(110) surface at half-filling and U = 8. SOPT-HF 
results for the first three layers from the surface. Dashed 
lines: local approximation for the self-energy. Solid lines: 
non- locality of self-energy fully taken into account. 



Compared with the free DOS in Fig. 1, the differ- 
ences between the QDOS for the different layers are 
much smaller in general. This is a consequence of quasi- 
particle damping due to the imaginary part of the self- 
energy. Generally this indicates a damping of the surface 
Friedel oscillations due to electron correlations. Note 
that the effect is also reproduced by the local approxi- 
mation. The strong periodic deformations of the DOS 
in the vicinity of the surface that could be observed for 
the free (U = 0) system (Fig. 1), are now confined to 
the resonance in small energy range around the Fermi 
energy. Here they appear as tiny wiggles only. 

The effects that result from the reduction of the co- 
ordination number at the surface show up more clearly 
when comparing results for different surfaces. For this 
purpose we also considered the (110) and the (111) sur- 
face of the sc lattice. Thinking the respective system 
to be built up from layers parallel to the surface, each 



site except for those in the first layer has 



(in) 



'II 



0) nearest neighbors within the same and 



,(no) 



4 (z 



(in) 



6) nearest neighbors in the adja- 



cent layers. This implies z^g 10 



4 and 



y (iu) 



3 for 



the coordination numbers of the top-layer sites, and we 
havezf 0) >4 110) >4 m) . 

This sequence can be recovered in the results for the 
top-layer self-energy for the different surfaces. Fig. 9 
shows their real parts for U — + and bulk band-filling 
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Fig. 12. The same as Fig. 11, but for the sc(lll) surface. 



n = 0.6. The comparison shows significant differences 
which prove the top-layer self-energy to be sensitively 
dependent on the type of the surface. The more open 
the surface, the stronger is the narrowing of the top- 
layer free DOS and thus the narrowing of the sclf-encrgy. 
In Fig. 9 this effect can be seen most clearly at high 
energies. 

It is even more apparent in Fig. 10 where the cor- 
responding imaginary parts are shown in comparison 
with the imaginary part of the bulk self-energy. The 
strongest differences with respect to the bulk are found 
for the (111) surface. As has been mentioned above, 
the lowered variance of the top-layer free DOS implies 
a reduced charge density at the surface: n s < n < 1. 
We find Us 1 = 0.38 and ni 110 "* = 0.52 to be compared 
with rig 100 ' = 0.56 and n — 0.6. According to Eq. ( |lg| ) 
this results in a diminished total weight of the imaginary 
part of the top-layer infinitesimal self-energy. Again the 
weight is reduced in the occupied part E < mainly. 

In the case of the (100) surface, differences between 
the local self-energy at the surface and in the bulk are 
worth mentioning for the top layer only (as discussed 
above). The same holds for the (110) surface. In the 
case of the (111) face there are non- negligible differences 
also for the second layer, while the local self-energy of 
the third (and all other) surface layer is almost indistin- 
guishable from the bulk self-energy. 

Finally let us investigate the effect of the non-local 
parts of the self-energy on the QDOS. The local QDOS 
in the first three layers for U = 8 and at half-filling 
(n = 1) is shown in Fig. 11 for the (110) and in Fig. 12 



for the (111) surface (solid lines). In both cases we find 
the spectra for all layers to consist of the two high-energy 
charge-excitation peaks as well as of the resonance at the 
Fermi level. Compared with the results for the (100) sur- 
face (Fig. 8), the fine structure within the resonance is 
much stronger. This has to be regarded as a pure sur- 
face effect since corresponding structures are found in 
the free DOS at the respective surface as well. The nar- 
rowing of the top-layer QDOS is another surface effect 
that is more pronounced for the (110) and even more for 
the (111) than for the (100) surface. This is consistent 
with Eq. plj ) when taking into account the coordina- 
tion numbers of the top-layer sites at the respective sur- 
faces. Due to the imaginary part of the self-energy the 
differences between the QDOS for different layers are 
diminished to a large extent compared with U = 0. If 
this can be extrapolated to the asymmetric case n ^ 1, 
at least qualitatively, it implies a damping of the surface 
Friedel oscillations in the layer-dependence of the charge 
density. 

The QDOS shown in Fig. 11 and Fig. 12 has been cal- 
culated including the non-local parts of the self-energy 
up the third neighbor shell (solid lines). This is com- 
pletely sufficient to obtain fully converged results. For 
s = 1, i. e. including the on-site and the nearest-neighbor 
non-local part only, we get almost the same results. It 
can thus be concluded that for all surfaces considered 
the non-local terms beyond nearest neighbors are quite 
unimportant as concerning their effects on the QDOS. 
Differences to the fully converged results show up in the 
case s = 0, i. e. for the local approximation (dashed 
lines). However, even for the more open (110) and (111) 
surfaces the local approximation qualitatively recovers 
all features in the spectra. Concerning the top-layer 
QDOS, it is slightly worse than for the (100) surface. 
For the second- and the third-layer QDOS the errors 
are comparable to those for the bulk case. Again we 
find an overestimation of the weight of the resonance 
around E = while the weight of the charge-excitation 
peaks is slightly too small. 

IV. CONCLUSIONS 

The behavior of the electronic self-energy and of the 
quasi-particle density of states at the low-index surfaces 
of a simple cubic lattice has been investigated within 
the framework of the semi-infinite Hubbard model. Our 
main interest has been focussed on the non-locality of 
the self-energy and its effects on the electronic struc- 
ture. Within the weak-coupling approach that has been 
employed for the present study, the non-locality can 
fully be taken into account. The results for infinitesi- 
mally small interaction U are exact. For finite U and 
at half-filling the SOPT-HF approach is expected to 
yield qualitatively reliable results at least. SOPT-HF 
studies have been performed previously for infinitely ex- 
tended lattices with full translational symmetry (e. g. 
Refs. 0,||). 

Let us briefly summarize the results found for the be- 



havior of the self-energy Y,\-'{E) at a single-crystal sur- 
face: The reduced translational symmetry implies that 
there are inequivalent sites i and j within the same 
neighbor shell. This results in a directional dependence 
of the non-local self-energy parts. The directional de- 
pendence has been found to be weak and to be confined 

(2) 

to the very surface. The local self-energy T, u (E) at 
the surface differs from its bulk value for the top layer 
only (in the case of the open (111) surface there are 
non- negligible differences also for the second layer) . Dif- 
ferences of the non-local parts (i =/= j) compared with 
the bulk case are significant if there is at least one site 
within the top layer. As a consequence of the reduced 
coordination number of sites at the very surface, there 
is a narrowing of the free top-layer DOS. This results 
in a corresponding narrowing of the local top-layer self- 
energy. Furthermore, the reduced coordination number 
results in a diminished (enhanced) surface charge den- 
sity below (above) half-filling which (in each case) im- 
plies a lowered total weight of the imaginary part of the 
top-layer self-energy. The decrease of the absolute mag- 
nitudes of the non-local parts with increasing distance 
between the sites i and j is as fast as for the bulk, but 
considerable faster than for the two-dimensional case. 
All effects mentioned have been found to be the stronger 
the more open the surface. 

The effects of the non-locality of the self-energy on 
the quasi-particle density of states can be summarized 
as follows: Qualitatively, all effects that show up in the 
layer-dependent QDOS at the surface are the same as 
those for the bulk QDOS: The inclusion of more and 
more non-local self-energy parts reduces the weight of 
the quasi-particle resonance at the Fermi level and pro- 
nounces the high-energy charge-excitation peaks. Al- 
though the top-layer self-energy is sensitively dependent 
on the type of the surface, the effects of its non-local 
parts on the QDOS arc of the same order of magnitude 
for the different surfaces considered. In all cases the lo- 
cal approximation is found to yield qualitatively correct 
results. The errors introduced by the local approxima- 
tion in the surface region are comparable to those in the 
bulk. With respect to the QDOS the local approxima- 
tion is only slightly better than for the two-dimensional 
case. However, there is one main difference that dis- 
tinguishes between the D — 2 lattice on the one hand 
and the semi-infinite (D = 3) lattice on the other: In 
the latter one gets nearly convergent results by taking 
into account merely the on-site and the nearest-neighbor 
self-energy parts. On the contrary, up to s « 5 shells 
have to be included to obtain convergence for D = 2. 

It depends on the type of the physical problem to be 
investigated whether or not one can tolerate the errors 
introduced by the local approximation. This stresses the 
need for a quantitative estimate of the effects of non- 
locality which has been given here for the semi-infinite 
Hubbard model. Our study should have shown that if 
one agrees on the use of the local approximation for the 
D = 3 bulk, one may likewise consider a surface without 



the necessity for further justification. In particular, this 
is important for the application of the dynamical mean- 
field theory (DMFT) §, § to the semi-infinite Hubbard 
model. This means to apply those many-body tech- 
niques that have been developed for the (approximate 
or even exact) solution of the infinite-dimensional Hub- 
bard model. For D = oo the local approximation be- 
comes exact, and there is a universal functional relation 
between the local self-energy and the QDOS. Within the 
context of the DMFT the main question thus concerns 
the effects of the non-local parts of the self-energy on 
the QDOS which could be estimated within the present 
study. 
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